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Abstract 

In this paper, we investigate the holographic phase transition with dark matter sector in the 
AdS black hole background away from the probe limit. We disclose the properties of phases mostly 
from the holographic topological entanglement entropy of the system. We find the entanglement 
entropy is a good probe to the critical temperature and the order of the phase transition in the 
general model. The behaviors of entanglement entropy at large strip size suggest that the area law 
still holds when including dark matter sector. We also conclude that the holographic topological 
entanglement entropy is useful in detecting the stability of the phase transitions. Furthermore, we 
derive the complete diagram of the effects of coupled parameters on the critical temperature through 
the entanglement entropy and analytical methods. 
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I. INTRODUCTION 


The AdS/CFT correspondence provides us a powerful approach to holographically study strongly interacting 
low energy physics in condensed matter systems. According to this correspondence, the d dimensional strongly 
interacting theories on the boundary are dual to the d+I dimensional weakly coupled gravity theories in the 
bulk iHSl- The most simple holographic superconductor model dual to gravity theories is constructed by 
applying a scalar field and a Maxwell field coupled in an AdS black hole background ^-l6|. Since then, a lot 
of more complete holographic superconductor models were also taken into account, such as the holographic 
superconductor models in Einstein-Gauss-Bonnet gravity, Horava-Lifshitz gravity, non-linear electrodynamics 
gravity and so on. These typical examples have attracted considerable interest for their potential applications 
to the condensed matter physics, see [7|-[3^. 

Recently, a gravity theory with dark matter sector was proposed in [3^, l3^ ■ This new gravity theory was 
also considered in holographic superconductor models, which was constructed with a scalar field, a Maxwell 
field and another additional U(l)-gauge field corresponding to the dark matter one 


37 


Liouville eigenvalue and matching semi-analytical methods, it has been disclosed in 


. With Sturm- 


37j| that the dark matter 


sector can bring rich physics in the new holographic model. Very surprisingly, this new model also allow 
superconducting solutions corresponding to retrograde condensation. In order to further study the effects of 
the dark matter sector on the holographic phase transition from other aspects of the scalar operator, the free 
energy and so on, we will have to turn to the numerical methods. 


The papers 


33 . |40j have shown us a novel way to calculate the holographic entanglement entropy of a 


strongly interacting system from a weakly coupled gravity dual according to the AdS/CFT correspondence. 
In this way, the holographic entanglement entropy has recently been applied to study the properties of phase 


transitions in various holographic models 




57j . The entanglement entropy representing the degrees of 


freedom of the systems turns out to be a good probe to investigate the critical temperature and the order 
of the holographic phase transition. It was argued that the discontinuous slops imply the second order 
phase transition and the jump of the holographic entanglement entropy corresponds to the first order phase 
transition. It is meaningful to examine whether the holographic entanglement entropy approach is still useful 
in the holographic superconductor model with dark matter sector. From the other aspect, it was found that 
the entanglement entropy and thermal entropy behave qualitatively the same for large width strip 53|, l54| . In 
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accordance with the area law, it was found that the holographic entanglement entropy goes linearly for large 
strip 53- Then it is expected that we can test the area law with dark matter sector from the entanglement 


entropy side. At last, the retrograde condensation phenomenon was observed in 


37l| , which usually corresponds 


to unstable solutions 60|, l 6 l| . We will also try to disclose the stability of the retrograde condensation solutions 
through the holographic entanglement entropy method. 

The next sections are organized as follows. In section II, we review the construction of the holographic 
superconductor model with dark matter sector in the four dimensional AdS black hole spacetime beyond the 
probe limit. In section III, we study the properties of the holographic phase transitions by examining in detail 
the behaviors of the holographic entanglement entropy. We also give some analytical understanding of the 
phase transition properties. We summarize our main results in the last section. 


II. EQUATIONS OF MOTION AND BOUNDARY CONDITIONS 


The holographic superconductor model with dark matter sector is constructed by a scalar field and two 
gauge fields coupled in the AdS black hole background. The generalized Lagrange density of 4-dimensional 
spacetime with dark matter sector reads 37| : 

C^R+^- ^[^F^^Fmn + IVmV' - ^ ^ ^b^^Fmn], (1) 

where tpir) is a complex scalar field with mass m. Am stands for the ordinary Maxwell field and Bm is the 
additional U{\) gauge field representing the dark matter one. —3/L^ is the negative cosmological constant, 
where L is the AdS radius which will be scaled unity in our calculation. 7 describes the backreaction of 
matter fields on the background. When 7 —> 0, we return to the holographic model in the probe limit 241. 
a is the coupling parameter between the two gauge fields, which should be small and on the order of 10 “^ 


according to present astronomical observation 


58l| . In fact, when studying holographic superconductors, one 


is normally interested in the CFT on the boundary. It is not an issue whether the dual bulk gravitational 
theory is realistic or not. So in this paper, we study a varying in large range a G [0, 2.5]. 

The Einstein equations for the system can be written in the form 


Rmn — -^gMNR — SgMN = ■^'yTMN, 
where Tmn is the energy-momentum tensor expressed as 

Tmn = Fm^F^ + Bm^B^ + uBm/sF^ + 2V Nip + “^AmAmiP"^ 

+gMN{—^FMNF^^ — -BmnB^^ — —FmnB^^ — VmV’VmV’ — AmA^iP"^ — m’^ip'^). 


( 2 ) 


(3) 
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With the variation of the matter fields, we get the corresponding equations of motion: 


= 0 , 


VmV“'0 — — m^ip = 0 


VmB^^ + = 0 . 


Putting (6) into (4), we get the equation 


Y7 TpMI\ 

Vm-t — 


MN ‘ 2 .pj‘^A^ 


where ce = 1 — 


Putting (7) into (6), we arrive at 


= 0 , 


In this work, we simply take the metric solutions and other matter fields in the forms: 


ds^ = -g(r)e-^FUf + ^+r^idx^ + dy^), 

9[r) 


A = (l){r)dt, B = ri{r)dt, ip = ip{r). 


Then the Hawking temperature of the black hole is expressed as 


where r+ is the horizon of the black hole satisfying g{r+) = 0. We also need x(r —>■ oo) = 0 to recover the 
AdS boundary. 

From above assumptions, we can obtain the equations of motion as: 


\r 2 g J g g^ 


\ r If ao 


- + ^ ]v + -—(p = 0 , 

r 2 J a g 


'- (2a -i. 

I L2 r 


+ K^r 2gip'‘^ + + 7 ?'^ + acp'r]') + 2rn?ip‘^ ^ 2e (p ip _ 

9 
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x' + 






= 0 , 


(16) 


Where = j. When a hairy black hole with V' ^ 0 appears, we have to solve these equations by numerical 
methods. Near the AdS boundary (r —>■ oo), the asymptotic behaviors of the solutions are 

■0- , V'+ 


0 ■ 


A I -, - 

rj^/\ — 


f X 0) 


(17) 


with A± = (3 ± -09 + 4 to^)/ 2, where fj. and p can be interpreted as the chemical potential and charge density 
in the dual theory respectively. The other two operators ^ and w are dual to the U(l) gauge field ??(r). 

At the horizon of the black hole, we can impose proper boundary conditions as: 


0(?') = 00 + 01 (’^ -r+)-\ -, 

(18) 

0(r) = (j)oir - r+) + 0i(r - r+f H-, 

(19) 

Vir) = Vo{r - r+) + pi{r - r+)^ H-, 

(20) 

9{r) = goir - r+) + gi{r - r+)'^ H-, 

(21) 

Xir) = Xo + Xiir - r+) -, 

(22) 


where the dots denote higher order terms. Putting these Taylor expansions into equations, we are left with 
five independent parameters r+, 0 oj 00 ; bo and xo at the horizon. The scaling symmetry 

r ^ ar, t ^ at, 0 — >■ a0, rj — >■ ap, g — >■ a^g, (23) 


can be used to set r+ = 1 . 

Choosing = — 2 > — | above the BF bound 


1 


the second mode 0+ is always normalizable. To get a 


stable theory, we will fix = 0 and use the operator =< 0+ > to describe the phase transition in the 
dual CFT. For different values of 0o, we can rely on the parameters 0o, bo and xo to search for the solutions 
with the boundary conditions 0_ = 0, x(^ oo) = 0 and 0 fixed as an constant. We mention that the 
solutions go back to the case without dark matter sector when ^ 0. In the following discussion, we will 
divide the phases into two cases: 1 + + ^ > 0 corresponds to the thermodynamically stable solutions and 

1 + a— + % < 0 is the thermodynamically unstable phase transition. 
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In the case oi ip = 0, we get the exact analytic solutions in normal phase, a Reissner-Nordstrom-AdS black 
hole, which is given by 


9 = — - — + ^ -5-X = 0, (p = pi -), V = wi -), 


L2 ^ J.2 

where M can be interpreted as the mass of the black hole. 


r+ r 


r+ r 


(24) 


III. HOLOGRAPHIC PHASE TRANSITIONS IN ADS BLACK HOLE BACKGROUND 

A. The stable phases with 1 + 


In this part, we focus on the holographic entanglement entropy(HEE) of the phase transition satisfying 


l + a-A-^ > 0. The authors in Refs. 


39 


4011 have provided a method to compute the entanglement entropy 


of conformal field theories (CFTs) from the gravity side. We consider a belt geometry A with a finite width 
along the x direction and infinitely extending in y direction as: 0 <i/<L, where I is defined 

as the size of region A, and L is a regulator which is set to infinity. Minimizing the area of hypersurface 7 ^ 


whose boundary is the same as the stripe A, we can deduce the entanglement entropy for A as 

.2 1 1 


52| 


S = 


dz- 




(25) 


with 


1 = f" - 1 - 

2 Je V(4 - zAz'^gizY 

where z* satisfies the condition ^|z, =0 with z = ^ and r = -1 is the UV cutoff. 




(26) 


FIG. 1: (Color online) The holographic entanglement entropy as a function of the temperature T for fj, — 1, tA = 0.1, 
= —2 and I = 2. The dashed lines in each panel corresponds to the entanglement entropy of a pure AdS black 
hole and the solid lines correspond to cases of superconducting solutions. In the left panel, we have choose ^ = 1 and 
various a: a = 0 with Tc = 0.02533 (Red); a = 0.25 with Tc = 0.02292 (Green); a = 0.5 with Tc = 0.02079 (Blue). 
The right panel correspond to a = 0.25 and various 7 : 7 = 0.9 with Tc = 0.02497 (Red); ^ = 1.0 with Tc = 0.02292 

(Green); ^ = 1.1 with Tc = 0.02091 (Blue). 


We present the holographic entanglement entropy as a function of the temperature T in Fig. 1 with /r = 1, 
= 0.1, = —2, - > 0 and I = 2. The left panel is the cases of u = 1, = 0.1, nA = —2, - = 1, / = 2 
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with various a and right panel shows the cases of /r = 1, = 0.1, m? = —2, a = 0.25, I = 2 with various 


^ > 0. It can be easily seen from the pictures that when the parameters fixed, the entanglement entropy 
decreases monotonously as we choose a smaller Harking temperature T. It means there is a reduction in the 
number of degrees of freedom due to the condensate generated in the phase transitions 54| . 

Decreasing the tem pera ture, we should determine the physical curve by always choosing the point of lowest 


entropy at a given T 


52j| . For each set of parameters, we find a threshold temperature Tc, below which the 


hairy black hole appears. We also obtain approximate formulas for the holographic entanglement entropy of 
the normal state (Sbh) and superconducting state (Ssc) hr the case of /r = 1 , = 0 . 1 , = — 2 , I = 2 , 


a = 0.5 and 4=1: 


Sbh ~ -0.309 + 0.557(T - T^) + 3.434(r - T^f, 
Ssc ~ -0.309 + 0.862(T - TJ + 17.273(r - 


(27) 

(28) 


These formulas suggest that Sbh\t=Tc = Ssc\t=t^, \t=t^ ^ \ T=Ta and the jump of the slop of the 

entanglement entropy signal that some kind of new degrees of freedom like the Cooper pair would emerge in 
the new phase. 

In order to show there is a phase transition, we calculate the free energy of the system with /x = 1, = 0.1, 

rn? = —2, a = 0.5 and ^ = 1 in Fig. 2. The free energy of superconducting state lies below the free 
energy of the normal state suggesting that the superconducting solutions are thermodynamically stable. More 

calculations show that there are thermodynamically stable solutions for all 1 + + ^ > 0. 
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FIG. 2: (Color online) The free energy in cases of ^ = 1, = —2, = 0.1, a = 0.5 and ^ = 1. The solid blue line 

shows the free energy of the superconducting phase. The dashed line corresponds to the free energy of normal phases. 
The solid blue point corresponds to the phase transition temperature T = 0.02079. 



As shown with the solid blue point in Fig. 2, we find a critical phase transition temperature T = 0.02079, 
which is equal to the threshold temperature Tc = 0.02079 obtained from the behaviors of holographic entan¬ 
glement entropy. That means the holographic entanglement entropy can be used to search for the critical 
phase transition temperature. With fitting methods, we obtain approximate formulas for the free energy of 








normal state (Fbh) and superconducting state (Fsc) that: 


Fbh ~ -0.122 - 1.417(T - T^) - 10.136(T - - 38.412(T - 

Fsc ~ -0.122 - 1.417(T - T^) - 23.573(r - - 190.175(r - T^f. 


(29) 

(30) 


The front formulas suggest that: Fbh\t=t, = Fsc\t=t,, ^^^|r=Tc 7^ 

^ \t=Tc- la other words, the curve representing the physical phases with the lowest free energy is smooth 
at the critical temperature Tc. We conclude that the jump of the slope of holographic entanglement entropy 
corresponds to second order phase transition in the general holographic superconductor model with dark 
matter. 




FIG. 3: (Color online) The holographic entanglement entropy as a function of the belt I with p = 1, = 0.1, = —2 

and T = 0.018. The solid lines correspond to cases of superconducting solutions. In the left panel, we have choose 
^ = 1 and various a from bottom to top: a = O(Red), a = 0.25(Green), a = 0.5(Blue). The right panel from bottom 

to top corresponds to a = 0.25 and various ji- ~ 0.9(Red), ^ = l.O(Green), ^ = l.l(Blue). 


We also show the behaviors of the holographic entanglement entropy S with respect to the strip width I 
at a fixed temperature T=0.018 below the phase transition temperature Tc in Fig. 3. The solid colour lines 
denote the holographic entanglement entropy for superconducting phase. We see that for each line, S increases 
monotonically from a negative value to a positive value as we increase 1. Choosing a larger strip size I, the 
entanglement entropy becomes sensitive to a and s-H panels, the curves go linearly with I for large I, 


which means the area law holds in holographic models with dark matter sector 


52|. 


By applying the holographic entanglement entropy method, we go on to disclose in detail behaviors of the 
critical temperature Tc- We show the effects of a on the critical temperature Tc in Fig. 4. In the left panel, 
when we increase the value of a in cases oi > 0 and ^ = 1, Tc decreases slowly, which means that the 
phase transition is more difficult to happen. In contrast, in the right panel, with ^ = —1 and > 0, we find 
that Tc increase as a function of a or the phase transition is more easier to happen. Tc is always an constant 
Tc = 0.05875 for different values of a with = 0 in both panels. 











9 




FIG. 4: (Color online) The behaviors of the critical temperature Tc with p = 1, = —2 and the lines from top to 

bottom correspond to = 0 (Red), = 0.1 (Green), = 0.2 (Blue). The left panel represents the case of Tc as a 
function of a with ^ = 1 . The right panel shows the case of Tc as a function of a with = —1. 




k- 

FIG. 5: (Color online) The behaviors of the critical temperature Tc with /r = 1, = —2. The lines in the left panel 

from top to bottom correspond to a = 0.25 and various = 0 (Red), = 0.1 (Green), = 0.2 (^Blue). The lines 

in the right panel from top to bottom represent the cases: (1) a = 2.5, ^ = —1 (Red); (2) a —>■ 2, |- = —1 (Green); 

(3) a = 0.25, ^ = 1 (Blue). 


Now we turn to study how ^ and will affect the critical temperature T^. The curves in the left panel 
of Fig. 5 represent Tc as a function of ^ with /x = 1, vn? = —2, a = 0.25 and k^= 0, 0.1 or 0.2. It can be 
easily seen from the left panel that Tc decreases as the ratio becomes larger with > 0. We conclude that 
larger ratio — > 0 makes it more difficult for the scalar field to condense when considering the matter fields’ 
backreaction on the background. Here, Tc is also an constant Tc = 0.05875 when neglecting backreaction. 

In the right panel of Fig. 5, the bottom solid blue line shows that heavier backreaction corresponds to 
a smaller Tc in the case of /i = 1, = —2, a = 0.25 and = 1. With more calculations, we conclude 

that larger makes the scalar fields more difficult to condense in all thermodynamically stable phases. We 
observe an constant Tc = 0.05875 in the case of /i = 1 , m? = — 2 , a —2 and = — 1 . As a complete analysis, 
we also show behaviors of the critical temperature of thermodynamically unstable phases in the case of /i = 1, 
= —2, a = 2.5 and - =—1 satisfying 1 +a-+ <0. It is surprising that a larger backreaction parameter 

fl fl fl 1—1 


leads to a larger Tc for unstable solutions. The former analytical results in 37| are included in the left panel 
of Fig. 4 and the right panel of Fig. 5. So we obtain richer physics through holographic entanglement entropy 

















10 


approach. 




FIG. 6: (Color online) The scalar operator as a function of T with ^ = 1, = 0.1, and m? = —2. The left panel is 

with ^ = 1 and various lines from bottom to top as: a = 0 (Red, Tc = 0.02533), a = 0.25 (Green, Tc = 0.02292) and 

a = 0.5 (Blue, = 0.02079). The right panel is with a = 0.25 and various lines from bottom to top as: ^ = 0.9 (Red, 

Tc = 0.02497), = 1.0 (Green, T, = 0.02292) and ^ = 1.1 (Blue, Tc = 0.02091). 


At last, we try to examine the correspondence between condensation gap and critical temperature 


17 |. 


We show the condensation of the scalar operator (0+)^/'^+ with = 1, = 0.1 and = —2 in Fig. 

6 . In the left panel, we exhibit the condensation of the scalar operator by choosing various a as: a = 0, 
a = 0.25, a = 0.5 and ■^- = 1. The increase of the parameter a develops deeper condensation gap. From 
another aspect, the increase of a as: a = 0, a = 0.25 and a = 0.5 corresponds to Tc = 0.02533, Tc = 0.02292 
and Tc = 0.02079 respectively. In the right panel, we choose different ratios of ^ and other parameters fixed 
to detect the behaviors of phase transitions. It is shown that larger ratios ^ develop deeper condensation 


gap and the different values of - = 0.9, = 1 and - = 1.1 correspond to the critical temperatures 

fl i-l (1 fl 


- = 0.9, -S- = 1 and - = 

MM ’ M M 

Tc = 0.02497, Tc = 0.02292 and Tc = 0.02091 respectively. With various parameters a and -, we arrive at the 
correspondence of a deeper condensation gap with a smaller Tc in the holographic model with dark matter 
sector. 


B. The unstable phases with l + + <0 


C pi 

In this part, we pay attention to the phase transitions satisfying 1 + + ^ < 0. We study the example 

p, = 1, = —2, = 0.1, a = 2.5 and ^ = —1 in Fig. 7. We show the holographic entanglement entropy 

in the left panel. Choosing the phases with lowest entanglement entropy, it is surprising that the condensed 
phase appears at high temperature T > Tc = 0.07231, which implies that this solution is unstable. This novel 


behavior was referred as retrograde condensation [37], l60j, . We find that there are solutions of retrograde 


p pi 

condensation phenomenon for all 1 + + ^ < 0. 


The free energy is powerful in studying the phase transitions. We calculate the free energy of the system 
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T T 


FIG. 7: (Color online) The phase transition in the cases of /r = 1, m? = —2, = 0.1, a = 2.5 and ^ = —1. The left 

panel shows the behavior of holographic entanglement entropy with I = 2. The middle panel represents the free energy 
of the system. In both panels, the solid red line corresponds to the superconducting phase and the dashed red line is 
with the normal phase. 


in the right panel of Fig. 7. It shows the free energy of this hairy black hole is larger than the free energy of 
the black hole in normal phase. Since the physical procedure corresponds to the phases with the lowest free 
energy, we arrive at an conclusion that the retrograde condensation superconductor solutions are thermody¬ 
namically unstable and the novel behaviors of the holographic entanglement entropy can be used to detect 
the thermodynamical stability of the phase transition. 


C. Analytical study of the condensation 


We would like to give an analytical understanding at the qualitative effects of a, ^ and on the phase 
transitions. At the phase transition point, we have tpir) = 0. Then we can rewrite the equations of motion 
around Tc in the form 



According to the behaviors of the equations, we can assume the solutions rjir) = -^(r), y = 0 and deduce 
the equations for (j)(r) and g{r) without dark matter term 




(33) 


9 


^ _9 
r 


+ ^r(j)'^ = 0 , 


(34) 


where 7 = K^(l-|-a^-|-^). We can taken 7 as the effective backreaction parameter to detect the changes of 
the critical temperature on the idea that a larger effective backreaction parameter 7 corresponds to a smaller 
critical temperature Tc 16|, [ 2 ^ . 
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When = 0, 7 is independent of a and This is in accordance with the fact that Tc keeps as an constant 
for different values of a and ^ in cases of = 0 in Fig. 4 and Fig. 5. If we choose a larger a with fixed ^ > 0 
and > 0, j increases and the critical temperature Tc becomes smaller. With — 1 < 0 and > 0, a 
larger a leads to a smaller 7 and a larger Tc- We have now covered all the qualitative properties of Fig. 4. 

In the left panel of Fig. 5, a larger ^ > 0 leads to a larger 7 and smaller Tc when a > 0 and > 0. 
Now we try to see how will affect the threshold temperature Tc. In the case of (1 + + jj^) > 0, a 

larger corresponds to a larger 7 and a smaller Tc. 7 is equal to zero for (1 + + It) = 0. In cases that 

(1 + a- + < 0, 7 becomes smaller and Tc becomes larger when we choose a larger Accordingly, the 

three lines in the right panel of Fig. 5 from bottom to top correspond to the decrease of (1 + + ^) as: 

2.25, 0 and -0.5 respectively. We argue that negative effective backreaction parameters with thermodynamical 
instability is a general properties in holographic phase transitions. 

IV. CONCLUSIONS 


We studied the behaviors of holographic metal/superconductor phase transitions in the presence of dark 
matter sector. We tried to explore the properties of the phase transitions by analyzing the holographic 
entanglement entropy of the system. It was showed that the entanglement entropy can be used to search the 
critical temperature and the jump of the slop of the holographic topological entanglement entropy corresponds 
to a second order phase transitions when including the dark matter sector. For larger strip size, the behaviors 
of holographic topological entanglement entropy suggest that the area law still holds with dark matter sector. 
We also found that the entanglement entropy serves as a good probe to the stability of the phase transitions. 
We obtained the thermodynamically stable conditions 1 -|- a- -I- % > 0 and our discussion are based on stable 

fl fl 

phases. In summary, we arrived at the conclusion that the holographic entanglement entropy can be used to 
explore the rich physics in holographic superconductor models with dark matter sector. 

We also derived the complete diagram of effects of the parameters a, and on the critical temperature 
Tc through the entanglement entropy methods. We found a and ^ don’t affect the critical temperature 
when neglecting the backreaction. The larger positive parameters ^ and make the phase transition more 
difficult to happen. When ^ > 0 and > 0, larger a makes the scalar fields harder to condense and 
larger a makes the scalar hair easier to form in the case of - < 0 and > 0. With various a and -, we 


also found that a smaller critical temperature corresponds to a 


obtained richer physics than the former analytical results in 


deeper condensation gap. In all, we have 


37|. By taking 7 = fc^(l -I- as the 
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effective backreaction parameter, the qualitative properties can be obtained through the analytical methods 
on the idea that a larger effective backreaction parameter 7 corresponds to a smaller critical temperature 
Tc- That means the field ri(r) and (^(r) coupled to determine the critical temperature. We argued that 
negative effective backreaction parameters developing thermodynamically unstable solutions may be a general 
properties in various holographic superconductor models. 
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